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ABSTRACT. Weighted voting games are frequently used in decision making. Each voter has a weight and a proposal is 
accepted if the weight sum of the supporting voters exceeds a quota. One line of research is the efficient computation of 
so-called power indices measuring the influence of a voter. We treat the inverse problem: Given an influence vector and 
a power index, determine a weighted voting game such that the distribution of influence among the voters is as close as 
possible to the given target value. We present exact algorithms and computational results for the Shapley-Shubik and the 
(normalized) Banzhaf power index. 

1. Introduction 

The European Economic Community, i. e. the predecessor of the European Union, was founded in 1957 by 
Belgium, France, Germany, Italy, Luxembourg and the Netherlands. Decisions were and are made via polling. 
More precisely, each country has a given voting weight, see Table and a proposal is approved if it is supported 
by at least 12, called the quota, out of the 17 votes. So Luxembourg possesses i G f a n votes. But let us consider 
the following thought experiment: The five countries except Luxembourg vote for their own without asking for 
Luxembourg's opinion. In any case the number s of supporting votes is even. Thus the vote of Luxembourg does 
not matter at all. For s > 12 the proposal will be approved and for s < 10 the proposal will be rejected. We may 
say that Luxembourg has absolutely no power in this weighted voting game. 



state 


votes 


population (2010) 


Germany 


4 


82144902 


France 


4 


62582650 


Italy 


4 


60017346 


the Netherlands 


2 


16503473 


Belgium 


2 


10783738 


Luxembourg 


1 


494153 



Table 1 . Voting weights and the quota in the European Economic Community of 1957. 



To have a measurement of the relative power of the members of voting schemes in more complex situations, 
power indices like the Shapley-Shubik or the Banzhaf power index were invented, see the next section for a defi- 
nition. Table|2]lists the values of the two mentioned power indices for three different weighted voting games. We 
remark that there is a non-linear dependency between weights, quota and the value of the considered power index. 

If we suppose that there exists some agreement on a fair distribution of power (<^i)i<j<n am °ri8 the n countries 
and a selection of an appropriate power index V, we obtain the inverse power index problem^] min \\d — V(x)\\> 
where \ is a weighted voting game for n voters and || • || an arbitrary suitable norm to measure the deviation, i. e. a 
rule to decide which assignment of voting weights is fairer. In practice the intended power distribution d may arise 
from the population numbers, either in the variant "One Person, One Vote", i. e. di is proportional to the population 
of country i, or by using Penrose square root law, i. e. di is proportional to the square root of the population of 
country i ||20ll22l . 



2000 Mathematics Subject Classification, Primary: 91B12; Secondary: 90B99. 

Key words and phrases, weighted voting game, simple game, complete simple game, Banzhaf index, Shapley-Shubik index, inverse 
problem. 

'Recent estimates for the population are taken from http://epp.eurostat.ec.europa.eu. 
^In II II the authors call this problem the voting game design problem. 
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state 


Wi 


SS 


BZ 


Wi 


SS 


BZ 


Wi 


55 


BZ 


Germany 


4 


0.233 


0.238 


4 


0.233 


0.229 


4 


0.233 


0.222 


France 


4 


0.233 


0.238 


4 


0.233 


0.229 


4 


0.233 


0.222 


Italy 


4 


0.233 


0.238 


4 


0.233 


0.229 


4 


0.233 


0.222 


Netherlands 


2 


0.150 


0.143 


3 


0.200 


0.188 


2 


0.100 


0.111 


Belgium 


2 


0.150 


0.143 


2 


0.050 


0.063 


2 


0.100 


0.111 


Luxembourg 


1 


0.000 


0.000 


1 


0.050 


0.063 


1 


0.100 


0.111 


quota 


12 






12 






11 







Table 2. Power indices for different weighted voting games. 



1.1. Related work. There is a vast literature how to compute the Shapley-Shubik index and other power indices 
in various circumstances, either exactly or approximatively. Few results are known about the inverse power index 
problem so far. Leech proposes a certain kind of fixed point algorithm and reports that it works considerably well 
in practice, whenever the number of voters n is not too small lfT6l [TBI , see also lfT4l . The author argues that for 
large n one may assume that the power index smoothly depends on the voting weights and thus Brouwer's fixed 
point theorem can be applied to prove the convergence of this approach. Another heuristic is described in (4|. 

Besides using the finiteness of the set of weighted voting games, the first general bounds stating that some power 
distributions can not be approximated by Banzhaf vectors too closely are given by Alon and Edelman [ 1 1. 

Since for each number of voters n there is only a finite set of weighted voting games one may simply solve the 
inverse power index problem by looping over the whole set. The enumeration of weighted voting games dates back 
to at least 1962 [18], where up to 6 voters are treated. For n — 7, 8 voters we refer e. g. to ll26l [T9l fl2ll . Bart 
de Keijzer presents a promising graded poset for weighted voting games in his master thesis ifTOl . see also ifTTl j^] 
n = 9 voters where successfully treated in l"P3l l23l . 

A similar inverse problem was solved in [3 1. 

1.2. Our contribution. We present the first exact algorithm, besides complete enumeration of weighted voting 
games, for the inverse power index problem. We computationally show that for n < 15 voters the Shapley- 
Shubik vector of every simple game, complete simple game or weighted voting game has a || • ||i-distance to 
(0.75, 0.25, . . . ) of at least |. For the Banzhaf power index the (conjectured) bound for the worst possible 
approximation in || • ||i-distance is a bit more involved. But it seems that for n > 2 voters (0.75, 0.25, . . . ) is 
one desired power distribution where the worst case is attained. We present computational results for n < 11 for 
simple games and n < 16 for complete simple games and weighted voting games. For more realistic instances we 
consider subsets of the European Union and determine the exact solutions for the inverse power index problem for 
small numbers of voters. Additionally we present some results on the set of achievable power distributions for the 
Banzhaf and the Shapley Shubik power index for n < 9 voters. 

1.3. Outline of the paper. In Section|2]we state the definitions of the three voting methods simple games, complete 
simple games and weighted voting games. We further introduce the Banzhaf and the Shapley Shubik power index. 
The known a priori estimates are the topic of Section [3] We give precise conjectures about the, with respect to 
approximation via Banzhaf and Shapley Shubik vectors, worst case power distributions. In Section [4] we give an 
exact algorithm for the inverse power index problem for simple games, complete simple games and weighted voting 
games for the Banzhaf and the Shapley Shubik power index via an ILP formulation. Another exact approach via 
exhaustive enumeration is described in Section [5] A promising combination of both methods as a branch&bound 
approach is outlined in Subsection |5.1| Computational results evaluating the proposed exact algorithms for the 
inverse power index problem and proving the stated conjectures for small numbers of voters are given in Section[6] 
To get a more global understanding we consider the sets of achievable Banzhaf and Shapley Shubik vectors in 
Section [7J and close with a conclusion in Section [7J Additionally we classify the complete set of desired power 
distributions which have the worst case approximation property for very small n in an appendix. 



"'We would like to remark that the counts for weighted voting games with 6 < n < 8 voters are wrongly stated in 1101 . but the methods 
should work. 
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2. Voting methods and power indices 

Consider a yes-no voting system for a set N = {1, 2, . . . , n} of n voters. The acceptance of a proposal depends 
on the subset of the supporting voters. In general this can be described by a Boolean function \ '■ % N — ► {0, 1}- A 
quite natural requirement for a voting system is monotonicity, i. e. x(U) < x(U') f° r a ^ U C U' C N. Monotone 
Boolean functions \ '■ 2 N ~ ► {0; 1} with xW = an d X(N) = 1 are called simple games. 

Definition 1. For a simple game x '■ % N {0, 1} a subset U C N is called a minimal winning coalition if 
x{U) = 1 and x{U\{i}) = 0/ or a ^ i ^ U. Analogously a subset U C N is a maximal losing coalition if 
X{U) = and x(U U {i}) = 1 for all ie N\U. 

Each simple game is completely characterized by either the complete list W of its minimal winning coalitions 
or the complete list L of its maximal losing coalitions. Both the set of minimal winning coalitions and the set of 
the maximal losing coalitions are antichains. For each simple game x there is a dual simple game x d defined via 
X d (U) = 1 — x(N\U). A voter i £ N is called a null voter if he is not a member of at least one minimal winning 
coalition. 

In many applications voting systems are restricted to smaller sub classes of simple games, e. g. complete 
simple games arising from Isbell's desirability relation [9|: We write i □ j for two voters i, j £ N iff we have 

X^{i} U U\{j}j > x(P) for all j £ U C A r \{i}. A pair (N, x) is called complete simple game if it is a simple 
game and the binary relation □ is a total preorder and we abbreviate i □ j, j □ i by i □ j. To factor out symmetries 
we assume 1 □ 2 □ • • • □ n and write coalitions U C N as characteristic vectors u. We call u a winning vector iff 
X(U) = 1, otherwise we call u a losing vector. 




Figure 1 . The Hasse diagram for < on N = {1, 2, 3}. 



To state an analogue of minimal winning coalitions for complete simple games we need another partial ordering: 

fc 

For two coalitions u = ( u\ ... u n ) , v = ( v± ... v n ) £ {0, 1}" we write u ^ v iff we have u i — 

i=i 

k 

Vi for all 1 < k < n and abbreviate u < v, u ^ v by u -< v, u ^ w, u ^ w by u txi v. With this we 

i=l 

have x( u ) x( w ) f° r a H u < v. By we denote the set of all shift-minimal winning vectors, i. e. winning 
vectors which are minimal with respect to -<, Similarly we denote by L the set of all shift-maximal losing vectors, 
i. e. losing vectors which are maximal with respect to ^. We remark that each complete simple game is uniquely 
characterized by either W or L, see e. g. In Figure [T| we depict the Hasse diagram for X on N = {1, 2, 3} 
and see that there are exactly 8 complete simple games (or antichains) for n = 3 voters; 7 with one shift-minimal 
winning vector and one with two shift-minimal winning vectors. 

An even more important sub class of complete simple games are weighted voting games given by non-negatives 
weights Wi for all i £ N and a quota q. The corresponding monotone Boolean function x is obtained by 

X{U) = 1 <=>• ^2wi>q. 

Algorithmically we can check whether a given simple game x is a weighted voting game by solving a suitable 
linear program, see Section [5] 

In Table [3] we state the known numbers of simple games, complete simple games and weighted voting games 
for n voters, see Q [T3] [23] [25] . To define two of the most important power indices for these three types of voting 
methods we need another definition: 

Definition 2. For a voter i £ N and a monotone Boolean function x, on i-swing is a set U C 7V\{i} such that 
X(U) = 0andx(UU {i}) = 1. 
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nl2345 6 7 8 9 

#sg 1 4 18 166 7579 7828352 2414682040996 5.613043 • 10 22 > 10 42 

#csgl3 8 25 117 1171 44313 16175188 284432730174 
#wvgl 3 8 25 117 1111 29373 2730164 989913344 

Table 3. Number of simple games, complete simple games and weighted voting games for n voters. 



With this we can express the Banzhaf vector (BZ(x, 1), . . . , BZ(x, n)) as 

BZ{x, i) = — -\{UC N\{i} I U is an i-swing}| , 

771 

where m is the total number of swings. Also the Shapley-Shubik vector (SS(x, 1), ■ • • , SS(x, n)) can be expressed 
using swings via 

SS ^ = ^/ E |C7|! • (n — |?7| — 1)!. 

U is an z-swing 

Since a set U is an i-swing in a simple game x if an d on ly if it is an i-swing in the dual game x d , both the 
Shapley-Shubik vector and the Banzhaf vector coincide for pairs of dual simple games. For both power indices the 
value for a null voter is zero. Removing a null voter from a simple game with n voters gives a simple game with 
77—1 voters, where the power indices of the remaining n — 1 voters coincide. 

3. A PRIORI ESTIMATES 

Since heuristics can only find good solutions but are not capable of providing lower bounds we are considering 
a priori estimates. The first fact coming to mind is the finite number of weighted voting games with n voters. An 
upper bound of 2" ~"+ 2 is given in ifTTl . Respecting symmetry we conclude that for n > 2 there are at most 
2 n 2 -n+2 . n j < 2« 2 . n n different Banzhaf- or Shapley-Shubik vectors in [0, 1]". Decomposing [0, 1]" into 2™ 2 • 77™ 
subcubes of side lenght results in at least one empty subcube using the pigeonhole principle. Thus there are 
points in [0, l] n such that the nearest Banzhaf- or Shapley-Shubik vector is at least 2 „+i ra apart. 

Having a closer look at the swings one can obtain a slightly sharper bound. For 77 voters the number of swings 
lies between 77 and ([^J + l) • see 0. Using the fact that the numbers r/i of swings for voter i all have 

the same parity we can conclude that either BZ(x, i) — BZ(x, j) or we have 

\BZ(x,i)-BZ(x,j)\> 



(LfJ + i)-( Lt j + i) 

for voters i and j. 

Both bounds tend to zero when 77 tends to infinity. Recently Alon and Edelman [ 1 1 found a lower bound that 
does not tend to zero: 

Theorem 1. Let n > k be positive integers, let e < -j^p^ be a positive real, and let xbe a simple game for 77 voters. 

n 

If ^2 -B(x> i) ^ e > tnen there exists a simple game "x! for k voters such that 

\\B( X )-BW)\\ 1 = '£\B(x,i)-BW,i)\+ J2 I^(X,0-0|< ( ^ + 1)£ +e. 

i=l i=k+l ^ ' 

n 

Alon and Edelman have chosen k = 2 and e = 0.01 to deduce that for a simple game with ^ B(x, i) < 0.01 

i=3 

the two-dimensional vector (B(x, 1), B(x, 2)) lies within || • ||i-distance smaller than i of one of the achievable 
Banzhaf vectors, i. e. (1, 0), (0, 1), or (h, i). Thus for 77 > 2 voters the Banzhaf vector of every simple game has 
an II • I1 1 -distance of at least min (0.5- Tg,0.02) = 0.02 to the power distribution (0.75, 0.25, 0, . . . ). 

We may improve this argument a bit by choosing k — 2 and e — yg to deduce that the minimal || • ||i-distance 
to the power distribution (0.75, 0.25, 0, . . . ) is at least |, which is best possible using the general theorem. 
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Conjecture 1. For n > 2 voters consider the following desired power distribution: d\ = 0.75, d2 = 0.25, 
and di = for 3 < i < n. For each simple game (complete simple game or weighted voting game) x we 

n n 

have the inequalities ||<SiS(x) — = £ |iSiS(x,*) — di\ > g and \\BZ(x) ~ d\\i — £ \^^(x^) ~ di\ > 

t=l i=l 

\i 1 ; ; i ; / 2fc m _i : m = mod 2, 

yJ— l w/zere fc x = 1, k m = < 1 _ 1 , „ for m > 2 and h = 2, Z 2 = 5, « m = 

'r^i ^ 8fc m _i — 1 : m = 1 mod 2 

f 2Z m _i + 3 : m = mod 2, , . „ 
S , r, i , r, for m > 3. 

8« m _i — 2 : m = 1 mod 2 



m 2 3 4 5 6 7 

^ 1 = Q.5 | = 0.4 H » 0.39474 g w 0.37975 §§ 0.37937 ^ w 0.37846 
Table 4. Values of j 21 - the conjectured lower bound for the Banzhaf index. 



Using r(n) = - 1 ' ) 2 — — a closed expression for these numbers j^, see also Table Q is given by k m = 

7 - 22m ' 2+T( i : )+ ^ r(m) and / m = 2 2 "- 3 +^) + n - 22 '"' 2+T( :, )+1 ' 23T(m)m foi- m > 2. Thus we have lim ^ = 
|| w 0.378378 > |. 

Conjecture [T] looks rather specific. A first hint why d\ = 0.75, di — 0.25 might be hard to approximate is the 
fact that this point maximizes the || • || i-distance to the achievable power vectors (1, 0), (0.5, 0.5), (0, 1), both for the 
Shapley-Shubik and the Banzhaf index. For some small n > 2 we compute the points or regions which maximize 
the || • ||i-distance to the achievable power vectors in Appendix |A| So far we have not seen a geometric structure 
behind the values for the Banzhaf index, but computing these numbers for simple games and weighted voting 
games with n < 16 suggests that the rather structured recursion formulas for k m and l m might be more than a pure 
coincidence. Having Theorem [T] in mind one might conjecture that power distributions with many zeros are hard 
to approximate whenever they are hard to approximate without the zero values. A first step towards understanding 
the inverse power index problem might be the results on achievable hierarchies in simple games, see J8). 

Additionally we conjecture that there is a similar result as Theorem[T]for the Shapley-Shubik index. In the next 
sections we will present some exact algorithms to solve the inverse power index problem and use them to prove 
Conjecture[T]for small n. 

n 

Conjecture 2. For each desired power distribution di with di > and Y] dj = 1, where n > 3, there exist two 

i=l 

n 

weighted voting games (complete simple games or simple games) \, x' w i tn \\BZ(x)—d\\i — ^2 W^{Xi *) — di\ < 

i=i 

j^landWSSW-dh = £ \SS(^,i)-di\ < |. 

4. An ILP formulation for the inverse power index problem 

Our first exact approach to solve the inverse power index problem is to utilize an ILP formulation using I LOG 
CP LEX 11 as an ILP solver. Therefore let n be the number of voters and d = (di, . . . , d n ) the desired power 

n 

distribution with di > and £ di = 1. To model the monotone Boolean function \ we use binary variables xjj 

8=1 

which equal x(U) for all subsets U C N. In order to express the values SS(x, i) and BZ(x, i) of the power indices 
for voter i we additionally introduce binary variables y^u which equal one iff U is an i-swing, see Equation |5]l 
With this we can state the Shapley-Shubik value for voter i as 

SS(t)= Pl = ^- Yl (\U\\-(n-\U\-l)l)> Vi , v 

UCN\{i} 

and the number of i-swings as Sj = £ Vi,U f° r all 1 < i < n. To measure the deviation \ P i — di \ for voter i 

UCN\{i} 

we introduce variables Si and obtain the following binary ILP formulation for the inverse Shapley-Shubik index 
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problem with || • || i-norm on simple games: 

(1) min 

(2) s.t 
(3) 

(4) ft = i- £ (|ET|!.( 

UCN\{i} 

(5) 
(6) 
(7) 
(8) 
(9) 
(10) 

(11) 
(12) 

Of course we may directly eliminate all variables yijj and pi from this formulation. The equations and inequal- 
ities (|6j, 0, (|8j, and |9]i enforce that x is a simple game. To restrict \ t0 complete simple games only small 
adjustments are necessary: We replace the subsets U C N by its characteristic vectors u £ {0, 1}" and rewrite 
Inequality (|6| to x u < x v for all u, v g {0, 1}" with u < v, or more economically to 

(13) a; u <a; Vw € {0,1}", u u : |H|i ^ 0, 

where — > i> denotes the set of elements arising from v by an right-shift of one of the rightmost Is in v. Examples 
are given by -+(1,0, 1,0) = {(1,0, 0,1), (0,1, 1,0)} and -> (1,0,1) = {(1,0,0), (0,1,1)}. 

To restrict \ to weighted voting games we additionally have to introduce weights Wi > and a quota q > 0. 
To be able to decide whether w i _t 1 or S w i < 1 lVi tne framework of ILPs we restrict ourself to weights Wi 

which differ by at least one whenever they are different. (We may simply use integer weights, which could result 
in harder problems for the ILP solver.) To interlink the xy with the Wi and q we use 

(14) q - (1 - xu) ■ M - ^ Wi < VE7 C N 
and 

(15) -xu-M + Y^ m <q-l W C N, 



i=l 




Pi < + 5j 


VI < t < n, 


Pi > c?,: - Si 


VI < t < n, 


~ 1 /" r 1 1 M i „, 

" - 1^1 - • 


vl \ ?. \ 


2/i,!7 = £[/u{i} — a;*/ 


VI < i < n,l7 c 7V\{i}, 


Xrr Xtt\ :j\ 


V0 ^ U C TAT, 7 e U, 


xn, = 




X N = 1 




art/ G {0,1} 


VJ7 C TV, 


I/i,tr G {0,1} 


VI < i < n,f/ c N\{i}, 


Si > o 


VI < i < 71, 


ft > o 


VI < i < n. 



where M is an suitably large constant fulfilling M — 1 > ^uij. (We may choose A/ = 4n (^x^ ) <n+1 ^ 2 , see ifTTl 



i=l 



Theorem 9.3.2.1].) Now let us argue why Inequalities ( p"4| ) and ( p"5j ) are necessary and sufficient. If a?t/ = 1 in ( fT~4] > 



and a;t/ = in 



15 1 then the inequalities are equivalent to w i > 1 an d Wi < q — 1, respectively. Otherwise 



(for xjj £ {0, 1}) both sets of inequalities are valid due to M — 1 > ^ w^. So we have given ILP formulations for 

the considered three types of voting methods when the Shapley-Shubik index and the || • ||i-norm is used. 

For the Banzhaf index we have to tackle the fact that its value for voter i is given by the following non-linear 



equality BZ(i) = — and we have to evaluate the deviation \BZ(i) — di\ = di 

£ si <=i E 

n 

To this end we set s — Sj and Si = |s, — di ■ s\. Instead of minimizing the sum of the Si we introduce 

i=l 

n 

Si < a ■ s, where a £ (0, 2] is a properly chosen constant. With this we can state: If the following binary ILP 

has a solution then the corresponding simple game x approximates the desired power distribution d with an error 
in the || • || i-norm of at most a. Otherwise no such approximation is possible. Thus we can minimize the deviation 
in the || • l^-norm by performing a bisection on a. (Since s lies between n and ([^J + l) • (i a j +1 ), each two 

unequal Banzhaf vectors differ by at least (^^) 2 , and we only need 0(n) such steps.) 
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(16) y] Si = s, 

i=l 



(17) 


> S t < a ■ s, 

i=l 




(18) 


Si < di ■ s + Si 


VI < i < n, 


(19) 


Si >di-s-Si 


VI < i < n, 


(20) 


Si = ^2 Vi,u 


yi<i<n, 




UCN\{i} 




(21) 


Vi,u = xuu{i} ~ xu 


VI < i < n,U C 


(22) 


%u > xu\{j} 


V0 7^ £7 C AT,j e [7, 


(23) 


x% = 




(24) 


X N = 1 




(25) 


Xu e {0,1} 


VC7 C iV, 


(26) 


Vi,u e {o, 1} 


VI < i < n,U C 


(27) 


Si >o 


VI < i < n, 


(28) 


s t > 


VI < i < n. 



Of course we can apply similar adjustments for this ILP formulation as we did for the one modeling the Shapley- 
Shubik index to restrict the voting method x> described by the xu, to complete simple games or weighted voting 
games. 

4.1. Simplification for Conjecture [1] In order to computationally prove Conjecture [T] for the Shapley-Shubik 
index and small n we remove some variables and constraints from the previous ILP formulation. Since for the 

n 

power distribution d = (0.75, 0.25, . . . ) we have ^ $i = d\ + c?2 + (1 — Pi — P2) we can delete all variables and 

»=i 

constraints for i > 3: 

min di + d 2 + 1 — Pi — P2 

s.t. pi <di + Si VI < i < 2, 

Pi > di — Si VI < i < 2, 

^ = k\' (\U\\-(n-\U\-iy)j -y hU Vl<*<2, 

UCN\{i} 

Vi,u = xuu{i} ~ xu VI < i < 2, U C A^\{i}, 

> i^-j V0 ^ (7 C /V, C7 n {1, 2} ^ 0, j g 17, 
x = 0, 
xn = 1, 

a;c/G{0,l} VJ7C/V, 
y 4 ,c/e{0,l} VI < t < 2, 17 C JV\{»}, 

Pi,P2 > 0. 

5. An exhaustive enumeration approach for the inverse power index problem 

The conceptionally easiest exact algorithm to solve the inverse power index problem is to loop over the finite 
set of all voting methods and to pick that one that minimizes the deviation between the desired and the achieved 
power distribution. Table [3] proves this approach to be infeasible for simple games and n > 8. For complete 
simple games there exists a very effective parameterization, see [2 |. But, to our knowledge, nobody has determined 
the number of complete simple games for n > 10 yet. On the one hand the fraction of complete simple games 
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which are weighted voting games tends to zero as the number of voters tends to infinity, but on the other hand no 
parameterization of weighted voting games is known. Several recursive constructions, mostly using super classes 
of weighted voting games, are presented in the literature ifTTl fl2l 171 [191 . Here we will describe the approach of 
|UJ|6) which exhaustively generates the whole class of complete simple games in the first run and afterwards checks 
each complete simple game whether it is weighted or not via a linear program, see also E4l[T3ll . One variant is 
given by 



(29) y^u lWi >q VueW, 



(30) 1 + v i w i ^1 Vw e 

i=l 

(31) Wi > Wi + i VI < i < n — 1, 

(32) w n ,g>0. 

For n = 9 voters it seems inconceivable to solve 284432730174 linear programs in order to determine all weighted 
voting games, nevertheless this is - more or less - what was done in 1.23 jj A more effective way is described in 
lfl3l and uses partial complete simple games. 

As mentioned before a complete simple game can be characterized by its set W of shift-minimal winning vectors. 
Now we want to exhaustively generate complete simple games for a given number of voters using an orderly 
generation approach, see |2T1 . i. e. we start with an empty set W and add shift-minimal winning vectors, which are 
decreasing with respect to the lexicographical ordering <i e J5 ste P ste P- ^ n tne corresponding generation tree for 
a node W all successors W fulfill W CW 1 and y <i ex x for all x G W, y G W'\W. Now let L and V be the set 
of all shift-maximal losing coalitions corresponding to W and W. For a given set W ^ with lexicographically 

smallest element w we define L := |u G B n | 3u G W : v -< u\ U ju G B™ | $u G W U {x G B™ | x < kx w} : 

v >r U |, where B = {0,1} and condense the shift-maximal vectors of L to a set L C L. With this we have L C L 

and L C L' for all successors of W\ Thus whenever the previous linear program where is replaced by W and L 
is replaced by L is infeasible, we can prune the whole search tree below node W. 

Let us have an example for n = 7 voters to illustrate the proposed algorithm. For the partial set W = 
{1110000, 1101010, 1001111} of shift-minimal winning vectors we can determine the partial set 

L = {1101000, 1011010, 1100110, 1101001, 0101111} 

of shift-maximal losing vectors. Here 1101000 is an element of L since we have 1101000 -< 1110000 G W. Since 
the linear program corresponding to W and L is infeasible we know that none of the complete simple games given 
by the set W U X of shift-maximal winning vectors is weighted, where 

X c {0101111,0101110,0101101,0011111,0011110,0011101}. 

So at least 7 complete games are shown to be non-weighted by solving one linear program. 

Using the concept of orderly generation and partial complete games drastically reduces computation times. 
However the number of weighted voting games for n = 10 voters remains unknown. So this approach is limited to 
rather small numbers of voters. 

5.1. A branch&bound approach for the inverse power index problem for weighted voting games. Since 
we will see in Section [6] that the LP relaxation of our ILP models from Section [4] is rather weak we propose a 
branch&bound algorithm based on the discrete structure of partial complete games. 

As in Section|5]we start with an empty set W of the shift-minimal winning vectors. Again we generate complete 
simple games via orderly generation and prune the subtree of W whenever W can not be extended to a weighted 
voting game. To further prune the search tree we additionally use the information of W and the partial set L of shift- 
maximal losing vectors to determine some of the a;^-values in the ILPs of Section[4] We need a computationally 
cheap lower bound on the deviation between the values of the given power index for all voting methods x below 
node W and the desired power distribution d. If the best known solution so far is better than this lower bound, 
we can prune the search tree. A quite natural candidate for such a lower bound is the LP relaxation of our ILP 



^Additionally some computational tricks and heuristics were used to reduce the number of linear programs. 

5 Here one may also read the coalition vectors as integers written in their binary expansion and use the ordinary ordering < of integers. 
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models from Section [4] To know good solutions right in the beginning of our branch&bound tree we may utilize 
any heuristics for the inverse power index problem. In this paper we consider Leech's algorithm. 

6. Computational results 



state 


pop. 


state 


pop. state 


pop- 




Austria 


8200 


Germany 


82500 


Netherlands 


16300 


Belgium 


10400 


Greece 


11100 


Poland 


38200 


Bulgaria 


7800 


Hungary 


10200 


Portugal 


10100 


Cyprus 


700 


Ireland 


4100 


Romania 


21700 


Czech Republic 


10500 


Italy 


58500 


Slovakia 


5400 


Denmark 


5400 


Latvia 


3400 


Slovenia 


2000 


Estonia 


1300 


Lithuania 


2300 


Spain 


43000 


Finland 


5200 


Luxembourg 


500 


Sweden 


9000 


France 


60000 


Malta 


400 


United Kingdom 


60600 



Table 5. Population (in thousands) of the 27 member states of the European Union. 







simple games 




complete simple games 


inst. 


II -1 




nodes 


time 


ll-lli 


nodes 


time 


EUi 


0.00000 


• io-° 


1 


0.01s 


0.00000 • 10"° 


1 


0.01s 


EU 2 


7.69740 


• io- 2 


1 


0.01s 


7.69740 • 10~ 2 


1 


0.01s 


EU 3 


7.13793 


• io- 2 


5 


0.05s 


7.13793- IO" 2 


4 


0.01s 


EU 4 


6.30740 


• io- 2 


28 


0.01s 


6.30740 ■ 10~ 2 


9 


0.02s 


EU 5 


6.90250 


• io- 2 


303 


0.09s 


6.90250 ■ 10~ 2 


26 


0.04s 


EU 6 


4.18923 


• io- 2 


18265 


5.91s 


5.40190 ■ 10~ 2 


35 


0.11s 


EU r 


2.39402 


• io- 2 


71483741 


15.88h 


3.75078 ■ 10~ 2 


54 


0.41s 


EU 8 










1.78178- IO" 2 


210 


3.08s 


EU 9 










6.89922 • 10~ 3 


1271 


95.91s 


EU 10 










3.65220 • 10~ 3 


19520 


101.21m 


EU n 










1.31318- 10" 3 


349125 


55.25h 



Table 6. Results for the EU using the Shapley-Shubik power index. 



As a realistic input instance for the inverse power index problem we use the 27 countries of the European 
Union with population numbers popi given in Table [5] Instance EU n consists of the n countries with the largest 
population, where the desired power is given by 



di = 



3=1 

i. e. we apply Penrose square root law. In Table [6] and Table [7] we state the results for simple games, complete 
simple games, and weighted voting games, respectively, using the Shapley-Shubik power index. 

The weighted voting games corresponding to Table |7]can be represented by their quota and their weights yield- 
ing: [1; 1], [2; 1, 1], [2; 1, 1, 1], [3; 1, 1, 1, 1], [4; 1, 1, 1, 1, 1], [14; 5, 5, 4, 4, 3, 3], [18; 9, 8, 7, 7, 6, 5, 4], 
[41; 13, 11, 11, 10, 9, 9, 7, 6], [92; 30, 27, 27, 26, 23, 22, 17, 16, 14], [109; 34, 29, 29, 28, 25, 23, 18, 16, 14, 13], and 
[339; 96, 82, 81, 80, 71, 65, 51, 44, 38, 36, 36]. 

For n < 5 the optimal solutions for simple games and complete simple games coincide. For larger n simple 
games admit better approximations but require a huge amount of cpu time. For n < 9 the optimal solutions for com- 
plete simple games and weighted voting games coincide. We would like to remark that for weighted voting games 



the Inequalities (14i and (15i would suffice if the y-variables are binaries but adding the valid Inequalities (13i 



drastically reduces the necessary cpu time. 

We would like to remark that our ILP formulations for the inverse power index problem have a large integrality 
gap so that the lower bound stays zero for a long time during the solution process. Being more precisely we state 
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inst. 



weighted voting games 

nodes time lower bound 



EUi 


0.00000 


10" 





1 


0.01s 


0.00000 


10' 


-0 


EU 2 


7.69740 


10" 


-2 


3 


0.01s 


7.69740 


10 


-2 


EU 3 


7.13793 


10" 


-2 


9 


0.01s 


7.13793 


10 


-2 


EU 4 


6.30740 


10" 


2 


5 


0.02s 


6.30740 


10 


-2 


EU 5 


6.90250 


io- 


2 


23 


0.04s 


6.86700 


10 


-3 


EU 6 


5.40190 


io- 


2 


29 


0.17s 


2.52257 


10 


-2 


EU r 


3.75078 


io- 


2 


44 


0.40s 


4.97897 


10 


-3 


EU 8 


1.78178 


io- 


2 


165 


2.86s 


2.20380 


10 


-3 


EU 9 


6.89924 


io- 


3 


1449 


49.65s 


7.19156 


10 


-4 


EU 10 


5.35622 


io- 


3 


27517 


28.68m 


1.12706 


10 


-3 


EUu 


3.62854 


10" 


-3 


676774 


69.16h 


1.14349 


10 


-4 



Table 7. Results for the EU using the Shapley-Shubik power index. 



n 


1 


2 


3 


4 




5 




6 




7 8 


LP-bound 














1 


IO" 6 


1 


IO" 6 


1.4 


10" 6 


n 


9 


10 


11 


12 




13 




14 




15 


LP-bound 


9 • 10~ 7 t 


MO" 7 1 


10~ 7 3 


io- 7 


2 


io- 7 


8 


io- 7 


6 


io- 7 



TABLE 8. Lower LP-bounds for the instances EU n for simple games, complete simple games, 
and weighted voting games. 



the lower LP-bounds for the instances EU n for 1 < n < 15 players. It turns out that these bounds coincide for 
simple games, complete simple, and weighted voting games, but there are differences in the necessary computation 
times. The simple games for n > 13 took 14 seconds, 1 minute, and 6 minutes, respectively. For complete simple 
games the computation times for n > 13 are 2 seconds, 8 seconds, and 4 minutes. Solving the linear programs 
for weighted voting games with n > 13 took substantially longer: 4 minutes, 3 minutes and 55 minutes. For the 
instances hard n , see below, we remark that for all considered voting methods and 2 < n < 16 voters the lower 
LP-bound was given by zero and that we have observed a similar behavior of the computation times. 

Since the given LP formulation is that weak, in terms of the integrality gap, we might use the general restrictions 
on the achievable power vectors from Section [3] Another possibility to obtain a lower bound is to drop the x- 
variables from our ILP formulations and to aggregate the yi .[/-variables to Zjim- variables: 

n 

min Si 



s.t. d t - Si <Pi< d t + Si 

nl ■ Pi = j\ ■ (n - j - 1)! • y hj 

3=0 



VI < i < n, 
VI < i < n, 



(33) 



Z, n-1 < 



n — 1 
j 

n — 1 

n-l 
2 



Zij e N 
Si > 



VI < i < n, V0 < j < 



n-2 



Vl<i<n, n=l (mod 2), 
n — 2" 

VI < i < n, V0 < j < 



Vl<i<n. 
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Although the resulting lower bound, see the sixth column in Table 7^ is valid even for simple games, it turns out 
to be very useful, at least for small n. If we drop Equation (33 i the ILP decomposes into n simpler ILPs which can 
be solved more directly, e. g. via enumeration. 





simple games 


complete simple games 


inst. 


II - Hi 


nodes 


time 


II • Hi 


nodes 


time 


hard2 


0.5000 


1 


0.01s 


0.5000 


1 


0.01s 


hard3 


0.3333 


1 


0.07s 


0.3333 


1 


0.07s 


hard4 


0.3333 


6 


0.05s 


0.3333 


7 


0.01s 


hard 5 


0.3333 


14 


0.02s 


0.3333 


10 


0.02s 


hard 6 


0.3333 


28 


0.09s 


0.3333 


22 


0.03s 


hard7 


0.3333 


31 


0.11s 


0.3333 


35 


0.05s 


hards 


0.3333 


33 


0.19s 


0.3333 


43 


0.12s 


hardg 


0.3333 


43 


0.73s 


0.3333 


54 


0.33s 


hardio 


0.3333 


69 


2.30s 


0.3333 


127 


0.78s 


hardn 


0.3333 


183 


21.31s 


0.3333 


126 


2.85s 


hardi2 


0.3333 


177 


170.56s 


0.3333 


207 


18.62s 


hardi3 


0.3333 


151 


761.24s 


0.3333 


221 


141.79s 


hardi4 


0.3333 


2354 


4.32h 


0.3333 


194 


12.52m 


hardis 


0.3333 


3309 


23.49h 


0.3333 


353 


1.22h 


hard 16 








0.3333 


829 


53.67h 



Table 9. Results for the desired power distribution (0.75,0.25,0,...) using the Shapley- 
Shubik power index. 



Following Conjecture [T] and [2] for n > 2 we consider the instances hard n corresponding to the desired power 
distribution d = (0.75, 0.25, 0, . . . ). In Tableland Table[l0 we list the results for simple games, complete simple 
games, and weighted voting games, respectively, using the Shapley-Shubik power index. 





weighted voting games 


inst. 


II -Hi 


nodes 


time 


hard2 


0.5000 


1 


0.01s 


hard3 


0.3333 


3 


0.02s 


hard4 


0.3333 


6 


0.02s 


hards 


0.3333 


13 


0.03s 


hard 6 


0.3333 


19 


0.05s 


hard7 


0.3333 


31 


0.08s 


hards 


0.3333 


36 


0.17s 


hardg 


0.3333 


60 


0.43s 


hardio 


0.3333 


101 


1.47s 


hardn 


0.3333 


153 


4.05s 


hardi2 


0.3333 


208 


41.09s 


hardi3 


0.3333 


143 


252.71s 


hardi4 


0.3333 


339 


51.14m 


hardis 


0.3333 


579 


6.73h 


hard 16 


0.3333 


1087 


36.80h 



Table 10. Results for the desired power distribution (0.75,0.25,0,...) using the Shapley- 
Shubik power index. 



In Table 11 and Table 12 we list the results for simple games, complete simple games, and weighted voting 



games, respectively, using the Banzhaf power index. 



We may also add inequalities to ensure that all numbers of swings for a given voter have the same parity. 
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simple games 




complete simple games 


inst. 


ll-lli 


nodes 


time 


Mli 


nodes 


time 


hard2 


0.5000000 


1 


0.01s 


0.5000000 


1 


0.01s 


hard 3 


0.4000000 


6 


0.01s 


0.4000000 


5 


0.07s 


hard4 


0.4000000 


11 


0.01s 


0.4000000 


9 


0.17s 


hards 


0.3947368 


36 


0.02s 


0.3947368 


13 


0.25s 


hard 6 


0.3947368 


46 


0.10s 


0.3947368 


21 


0.26s 


hard7 


0.3797468 


177 


0.27s 


0.3797468 


15 


0.18s 


hard 8 


0.3797468 


1923 


7.80s 


0.3797468 


55 


1.15s 


hardg 


0.3793651 


11801 


82.46s 


0.3793651 


61 


2.42s 


hardio 


0.3793651 


167868 


57.60s 


0.3793651 


61 


17.11s 


hardu 


0.3784640 


311199 


140.85m 


0.3784640 


110 


131.67s 


hard 12 








0.3784640 


119 


18.07m 


hardi3 








0.3784399 


119 


87.85m 


hardi4 








0.3784399 


191 


9.48h 


hardis 








0.3783563 


217 


23.78h 


hardi6 








0.3783563 


248 


6.84d 



TABLE 11. Results for the desired power distribution (0.75,0.25,0,...) using the Banzhaf 
power index. 



inst. 



nodes 



weighted voting games 
time representation 



hard2 
hard3 
hard4 
hard 5 
harde 
hard7 
hardg 
hardg 
hard 10 
hardu 
hardi2 
hardi3 
hardi4 
hard 15 
hard 16 



0.5000000 1 0.01s [1;1,0] 

0.4000000 6 0.17s [3;2,1,1] 

0.4000000 26 0.16s [3;2,1,1,0] 

0.3947368 18 0.17s [5:4,1,1,1,1] 

0.3947368 45 0.28s [5:4,1,1,1,1,0] 

0.3797468 38 0.62s [14:11,3,3,3,2,1,1] 

0.3797468 55 1.41s [14; 1 1,3,3,3,2,1, 1,0] 

0.3793651 50 4.00s [33;26,7,7,7,5,2,2,1,1] 

0.3793651 73 24.20s [33:26,7,7,7,5,2,2,1,1,0] 

0.3784640 138 140.95s [80;63, 17, 17, 17, 12,5,5,2,2, 1,1] 

0.3784640 108 26.64m [80;63,17,17,17,12,5,5,2,2,1,1,0] 

0.3784399 111 75.33m [193;152,41,41,41,29,12,12,5,5,2,2,1,1] 

0.3784399 153 12.30h [193;152,41,41,41,29,12,12,5,5,2,2,1,1,0] 

0.3783563 151 1.31d [466:367,99,99,99,70,29,29,12,12,5,5,2,2,1,1] 

0.3783563 116 2.34d [466;367,99,99,99,70,29,29, 12, 12,5,5,2,2, 1,1,0] 



TABLE 12. Results for the desired power distribution (0.75,0.25,0,...) using the Banzhaf 
power index. 



7. Power distributions which are hard to approximate 

In the previous section we have computationally shown that the desired power distribution (0.75, 0.25, 0, . . . ) is 
hard to approximate by Shapley-Shubik or Banzhaf power vectors. Since it is reported that the proposed heuristics 
for the inverse power index problem perform well on practical instances one may conjecture that achievable power 
vectors are dense in those regions which are of practical importance and that the bad approximation property of 
(0.75, 0.25, 0, . . . ) is a very singular phenomenon. 

To shed some light on this matter we will consider two dimensional projections of the set of all achievable power 
distributions for small n. For 1 < i < j < n and a (normalized) power index V we define Afj = 

{(x, y) | 3 wvg x consisting of n ordered voters with V(i) — x,V(j) = y} . 

to be the set of all achievable power vectors for ordered weighted voting games consisting of n voters and a given 
power index V. By ordered we mean that voter i is as least as powerful as voter i + 1 for all 1 < i < n. 
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iff ' 




^ A 




^^^^^^^ '^^^fl 





FIGURE 2. Projections to the first two voters of achievable power distributions for 6 < n < 9 
voters and the power indices Banzhaf and Shapley Shubik, i. e. the sets A^ 2 f° r "P £ SS}. 



Lemma 1. For (x, y) E Afj with 1 < i < n we have 

1 > x > y > 0, 

(j — 1) • x + (n — j + 1) • y > l/o?" i = 1, anii 
i • x + (j - i) ■ y < 1. 

Proof. Since we assume "P(l) > 7^(2) > • • • > V(n) we have x > y and > V(j) for all ft, < j. Inserting 

n 

this into P(h) — 1 yields the stated inequalities. □ 

h=l 

In Figure [2] we illustrate these sets for 6 < n < 9 and T 3 G {SZ,55} in the following way. We have 
^= [0; 1] x [0, ij] and color the pixels of [0, 1] x [0, ^] in different colors. For (a;, y) is ruled out by Lemma T| 
we use red pixels and otherwise a scale of colors from dark blue to green. Here we use dark blue if (x, y) (actually 
all points in the near of (x, y)) is not contained in A^ 2 . By light blue we denote that some elements of A~i 2 aie 
in the near of (x, y) and by green we denote the case when there are many elements of A\ 2 WQ m tne near °f 
(x,y). The value for voter 1 increases from left to right and the value for voter 2 increases from top to bottom. The 
graphics on the left hand side of Figure [2] correspond to the Banzhaf power index and those on the right hand side 
to the Shapley Shubik power index. 

By a look at Figure |2]it seems that there might exist a certain subset of [0, l] 2 which remains empty (or dark 
blue) if n tends to infinity and we factor out the intrinsic discrete nature of the problem. But of course proving such 
a property might be a very challenging task. 

8. Conclusion 

With the inverse power index problem we have presented a very challenging optimization problem of practical 
importance. Besides exhaustive enumeration we have developed the first exact algorithm to tackle this problem 
via an ILP formulation. The first computational results are very promising but of course a lot of work using more 
sophisticated techniques from integer linear (or non-linear) programming has to be done. We give a very precise 
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conjecture on the power distributions with the worst case approximation property and reveal some structure of the 
set of achievable power distributions by the Banzhaf and the Shapley Shubik power index. In our opinion it would 
be very challenging to prove these conjectures. Another open problem, which might be not too hard to solve, is to 
prove an analogue of Alon's and Edelman's result for the Shapley Shubik or other power indices. We leave open 
the problem to design a practical algorithm for the inverse power index problem, i. e. an algorithm combining fast 
primal heuristics and techniques for lower bounds in order to solve the problem either exactly or with a certain error 
guarantee. 
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Appendix A. Regions which are hard to approximate 

In Conjecture [T] and Conjecture [2] we state that the desired power distribution d = (0.75, 0.25, 0, . . . ) is the 
worst possible example with respect to approximation in the || ■ ||i-norm using one of the three presented voting 
methods and the Shapley-Shubik power index. In this section we want to show that there are further desired power 
distributions which are equally hard to approximate, at least for small n. 

We restrict ourselves on complete simple games and weighted voting games, which coincide for n < 5 voters. 
For n < 2 and the Shapley-Shubik or the Banzhaf power index we only have the following achievable power 
distributions s\ = (1), s\ = (1,0), and s\ = (0.5,0.5) up to symmetry, i. e. we assume that the power is 
decreasing monotonely. Let d — (c?i, . . . ,d n ) be a theoretically possible desired power distribution, i. e. we have 

n 

di > and d^ = 1. To avoid symmetric situations we assume di > d i+1 for all 1 < i < n in the following. 

Since there exists only one such theoretically possible desired power distribution we cannot say something 
interesting for n — \ voters. 
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Lemma 2. Forn = 2 and the Shapley-Shubik or the Banzhaf power index the regions where min 1 1 d ~ V (x) 1 1 1 = h 

x 

are given by d = (0.75, 0.25). 

Proof. If d x > 0.75 we have |1 - di| + |0 - d 2 \ < 0.25 + 0.25 = \ and for \ < di < 0.75 we have d 2 > 0.25 
and|0.5-di| + |0.5-d 2 | < 0.25 + 0.25 = §. " " □ 

For n = 3 voters the achievable power distributions for the Shapley-Shubik index up to symmetry are = 

(l,0,0), S l = (| ! i,i), S l = (i,i,0),and S l=(i,i,i). 

Lemma 3. For n = 3 the regions where min \ \d — SS(x)\\i = | are given by 

• d! = f, i < d 2 < i 

• | < di < g, di + d 2 = 1 (i.e. d 3 = 0), 

• | < dx < §, d 2 = |, 

• di = |, j < d2 < |, and 

• ^ < di < i di + d 2 = | ("i.e. d 3 = \). 

Proof. (1) For di > | we have ||sf — d|L = 2 — 2di < |, where equality is attained iff di = |. 



(2.1) For | < d\ < | and d 2 < g we have [Jsj — d|| = 2d\ — | < |, where equality is attained iff di = |. 

(2.2) For | < di < | and d 2 > g we have d3 < g and ||s| — d|L = | — 2d3 < ^, where equality is attained 
iff d 3 = 0, i. e. di + d 2 = 1. 

(3.1) For I < di < I and d3 > g we have ||s 2 — d|L = | — 2di < |, where equality is attained iff di = |. 

(3.2) For 5 < di < |, d3 < g, and d 2 < | we have d 2 > | and ||s 2 — d|L = 2d 2 — | < |, where equality is 
attained iff d 2 = | . 

(3.3) For I < di < | and d 2 > | we have d 2 < | and ||s| — d|L = 1 — 2d 2 < |, where equality is attained 
iffd 2 = |. 

(4.1) For J2 < d\ < \ and d3 < | we have d 2 < \ and 1 1 sj| — d\\ = 2d$ < |, where equality is attained iff 
d 3 = |. 

(4.2) For ^ < d x < | and d 2 < | we have ||s| — d|| x = 2di — | < A, where equality is attained iff d x = ~. 

(4.3) For ^ < di < |, d 3 > |, and d 2 > | we have d 3 < | < | and ||s| - d]^ = | - 2d 3 < |, where 
equality is attained iff d3 = | . 

(5.1) For | < di < ^ and d 2 > | we have g < d3 < | and ||s| — d|| = | — 2d3 < |, where equality is 
attained iff d3 = | . 

(5.2) For I < di < ^ and d 2 < § we have ||s^ - d|| x = 2di - § < §. 

It remains to check, whether ||sf — d|L = | for the stated subspaces. This is easily verified. □ 

For the Banzhaf power index and 3 voters the achievable power distributions up to symmetry are b\ — (1, 0, 0), 

^=(!.|.|)^I = (|.i0),and6t=(|,i ) |). 

Lemma 4. For n = 3 the regions where min \\d — BZ(x)\\\ = § are given by d\ — |, ^ < d 2 < g ana" di = |, 
d 2 = |. 

Proof. (1) For di > | we have — d|| = 2 — 2di < |, where equality is attained iff d\ = |. 

(2.1) For I < di < I and d 2 > | we have d3 < | and — d|L = | — 2d3 < |, where equality is attained 
iff d 3 = 0. So here we can proceed as in the proof of Lemma |2j to deduce min \\d— BZ(x)\\ 1 < | for 

di^f. 

(2.2) For § < di < I and d 2 < g we have 1 1 fo^ — d|L = 2di — | < |, where equality is attained iff d\ = |. 

(3.1) For I < di < I and d 2 > ^ we have d 2 < | and ||6g — d| L = 1 — 2d 2 < |, where equality is attained 

iff d 2 = ^. For i < di < I and d 2 = ^ we have d 3 < | and — d|L = 4 < |. 

10 2 — ^ — 5 — 5 11 * 111 5 5 

(3.2) For i < dx < |, d 2 < ^, and d 3 < | we have d 2 > | and ||6| - dj^ = 2d 2 - | < |. 

(3.3) For \ < di < | and d 3 > | we have ||6| - d\\ % = | - 2d x < |. 

(4.1) For di < I and d 3 < | we have d 2 < | and ||6g — d|| = 2d3 < |, where equality is attained iff d 3 = |. 

(4.2) For I < di < ^ and d3 > | we have d 2 > | and — d|| = | — 2di < |, where equality is attained 
iffdi = §. 
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(5.1) For d x < I and d 2 > | we have di > |, A > d 3 > | and ||6f - d^ = | - 2d 3 < §. 

(5.2) For di < | and d 2 < § we have di > § and - d|| x = 2d x - | < §. 

For d 3 = g, see (4.1), we can assume d\ = d 2 = § due to (4.2) and can apply (5.1). It remains to check, 
whether — d|L = § for the stated subspaces. This is easily verified. □ 

For n — 4 voters the achievable power distributions for the Shapley-Shubik index up to symmetry are s\ = 

v ' ' ' /' *2 — (4' 12' 12' 12/' *3 — I3' 6' 6' u /' *4 — V 12 ' 4' 12' 12/' *S — l 2 ' 2' u ' u /' *6 — I2' 6' 6' 6/' 
„4 _ fJL JL JL M „4 _ ( 5_ 1 1 M ,4 _ /l 1 1 n \ „4 _ fl 1 1 1\ anr | „4 _ /l 1 1 1\ 
*7 — V 12' 12 ' 12' 12/' *8 — 1 12' 4' 4' 12/' *9 — I3' 3' 3' u / ' °10 — I 3 ' 3 ' 6 ' 6/' ^ll - U'4'4'4/' 

For the Banzhaf power index and 4 voters the achievable power distributions up to symmetry are fe 4 = (1, 0, 0, 0), 

(,4 _ fl 1 1 1\ Jk4 _ (3 1 1 n\ lA = (1 1 n q\ 1,4 _ H 1 I I\ ^4 _ (1 1 1 1\ ft 4 _ 
°2 - UO' 10' 10' 10 J' °3 - U' 5' 5' U J' °4 ~ b' 2' U ' U /' °5 ~ I 2 ' 10 ' 10 ' 10 / ' °6 - \ 2 i 6' 6' 6 J' °7 ~ 

f 5 1 1 1\ i4 _ f2 1 1 1 \ lA / 3 3 1 n i4 _ / 1 1 1 n) T, 4 — ( 1 1 1 1 "\ and /i 4 — ( 1 1 1 1 ) 

1 12' 4' 4' 12J'°8 - U' 5' 5' 5 / ' °9 ~ U' 8' 8' 8/'°10 ~ U' 3' 3' U /'°H _ V3' 3' 6' 6/' aIK10 12 ~ U' 4' 4' 4 J' 

Lemma 5. For n = 4 f/ze regions where min ||d — /B-Z^x)^ = | are g/ven fey di = g, d 2 = g and d\ = |, 

d 2 = |. 

Proof. (1) For d 4 > | we have ||6 4 — d|| = 2 — 2di < |, where equality is attained iff d\ = |. If d\ = | 
and c?2 < tl) then we have d 4 < d$ < ^ and ||6| — d]^ = g < |. If we have di = | and d 2 > jq then 
we have d2 < g, d± < d^ < and | &| — d L = 2d2 < §, where equality is attained iff d 2 = g and 
d3 = d4 = 0. 

(2.1) For | < di < | and d2 > g we have ||fe 4 — d|L = | — 2d3 < |, where equality is attained iff 
d3 = d 4 = 0. So here we can proceed as in the proof of Lemmaplto deduce min || d — BZ(x)\\ i < | for 

(2.2) For | < di < |, d 2 < g, and di + d± < | we have ||6g — d|L = 2di + 2d 4 — | < |, where equality is 
attained iff di + d 4 = | . 

(2.3) For | < di < |, d 2 < \, and di + d 4 > | we have d 4 < ^, d 4 > ^, and ||&| - d\\ ± < |, where 
equality is only possible for d 4 = so that (2.2) yields d\ — | and we can apply (1). 

(3.1) For | < di < | and d2 > we have d 2 < | and — d|| = 1 — 2d 2 < |, where equality is attained 
iff d 2 = j|. For 1 < dx < | and d 2 = ^ we have d 3 + d 4 < |, d 4 < i, and ||6 4 - d\\ l = | + 2d 4 < | 
where equality is attained iff d3 = d 4 = 4^ and d\ = |. But in the latter case we have — d = 0. 

(3.2) For \ < d 1 < |, | < d 2 < ^, and d 3 < | we have d 3 + d 4 = 1 - d 4 - d 2 > ^ - d 4 so that d 3 > ^ - ^ 
holds. Thus we have — d|L = | — 2di — 2d3 < ^ — di < |, where equality is attained iff d\ = ^, 
d2 = jy, d^ = d 4 = 4j, which equals 6|. 

(3.3) For | < di < | and d 3 > \ we have d 2 > | and - d|| x = | - 2di < §. 

(3.4) For | < di < | and d 2 < g we have d 3 < g and ||& 4 - d|| x = 2d 4 < | since d 4 < < g. 

(4.1) For d\ < \ and d3 + d 4 < g we have d 2 <\ and 1 1 b\ — d| 1 1 = 2d3 + 2d 4 < |, where equality is attained 
iff d3 + d 4 = |. For the latter case see (4.7). 

(4.2) For | < d 4 < | and d3 > | we have d2 > | and ||&| — d|L = | — 2d 4 < |, where equality is attained 
iff c?! = |. For the latter case see (5.1)-(5.4). 

(4.3) For di < |, d 2 < |g, and d 3 < | we have d 2 > |, d 4 < \, and ||6 4 - d\\ x = 2d 2 - | < |, where 
equality is attained iff d 2 = |g. In the latter case we consider x = 6| — d L. For d 4 < d3 < 4j we 
have x — ^ < | and for d3 > d 4 > we have a; = 2d3 + 2d 4 — A < |, where equality is attained 
iff d 3 + d 4 = |, i. e. d3 = d 4 = |, and 6^ — d|| < |. Finally, for d3 > 4^ and d 4 < 4^ we have 
x = 2d 3 - i < f . 

(4.4) For | < di < i and d 4 > g we have d 2 > d 3 > 1 and ||6 4 - d^ = 1 - 2di < |. 

(4.5) For | < di < \, d 2 > g, d 3 < |, and d 3 + d 4 > | we have d 4 <\ and ||bf - d|| x = | - 2d 3 - 2d 4 < |, 
where equality is attained iff d3 + d 4 = g. For the latter case see (4.7). 

(4.6) For | < di < \ and d 2 < g we have d 4 < d 3 < g and ||&| - d|| 1 = 2d 4 - | < |. 

(4.7) For | < d x < \ and d 3 + d 4 = | we have d 2 > ^, d 3 > 4^, d 4 < d 2 + d 3 < 1 - d x < | and 
||fo 4 — <2|| = 2d 2 + 2d 3 — | < |, where equality is attained iff d 2 + d 3 = |, i. e. d\ = |, d 2 = |, d 3 = g, 
and d 4 = 0. In the latter case we have ||6f — d|| < |. 
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(5.1) For d\ < |, d 2 > §, and d 3 < ^ we have d\ > |, d 2 < §, d 4 < g and i>| - d = \d\ — | + 

|d 2 - || + |d 3 - 1| + - || < ^ + M + h + I < I- 

(5.2) For d x < |, d 2 > |, and d 3 > ^ we have di > |, d 3 < |, and ||&f - d|| x = | - 2d 3 < |. 

(5.3) Fordx < §, d 2 < |,andd 4 < j| wehavedi > |§. For^ > | we have |[&f -d||j = 2d!+2d 4 -§ < § 
and for di < | we have ||&f — d|| = 2d 4 < |. 

(5.4) For ^ < di < |, d 2 < |, and d 4 > ^| we have d 2 > | and consider a; = — d . If d 3 > d 4 > | 
then we have x = | — 2d 2 < § for di > | and x = | — 2di — 2d 2 < §• If di < d 3 < \ then 
we have x = 2d\ — | < | for di > | and a; = for di < |. If d 3 > | and d 4 < | then we have 
x = 2di + 2d 3 - 1 < | for di > i and x = 2d 3 - | for di < |. 

(6) For j < di < ^ we have d 3 < d 2 < d 4 = 1 — di — d 2 — d 3 > l| and we consider x — ||6| 2 — d|| 
For d 4 < d 3 < d 2 < j we have a; = 2di — | < § , for d 2 > j and d 4 < d 3 < \ we have x = 
2di + 2d 2 - 1 < |, and for d 2 > d 3 > \ and d 4 < | we have a; = \ — 2d 4 < |. 
Finally we have to show that min ||d — BZ{x)\\ 1 = § for the stated regions. For d\ = 0.75, d 2 = 0.25 we have 

done this already by solving an ILP. Now let d\ = | and d 2 = |. We can easily check ||d — (6f ) || >= | for 
all 1 < i < 12 and \\d - BZ (b\)\\ 1 = §. □ 

We would like to remark that the weighted voting game x = 

[5; 4,1, 1,1,1] has (g, i, i, ^) as its 
Banzhaf vector and fulfills || (0.8, 0.2, 0, 0, 0) - BZ(x)\\ x < §. So the desired power distribution (0.8, 0.2, 0, 0, 0) 
is not an example for worst possible approximation property. So maybe for n > 5 the desired power distribution 
(0.75, 0.25, 0, 0, 0) is the unique example with worst possible approximation property for the Banzhaf power index. 
We remark || (0.75, 0.25, 0, 0, 0) - BZ{ X ) \\ x = §§ « 0.3947368. 

As far as the Shapley Shubik power index is concerned we have no idea if (0.75, 0.25, 0, . . . ) is the unique 
example with worst possible approximation property for large n or if there are regions which are hard to approxi- 
mate. In principle we can utilize our ILP formulation from SectionHto check whether a whole region like d\ = |, 
j2 — < \ can not be approximated with an error less than |l by letting the di variables with certain range 
bounds instead of numerical constants. As an example we have tried this for n — 6 and obtained the weighted 
voting game x = [5; 5, 1, 1, 1, 1, 1] approximating the desired power distribution d = (|, ^ , 0, 0, 0) with an 
erroTof\\d-SS(x)\\x = l<l. 
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